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1. INTRODUCTION 


The radio frequency bands below 10 GHz have been used 
for satellite-earth communications since the beginning of the 
space communications era. These bands are approaching saturation 
due to the ever increasing traffic. Consequently, for a number 
of years satellite communications engineers have been studying 
higher frequency bands for future satellite communications sys- 
tems to obtain an in-depth understanding of the propagation char- 
acteristics of earth-satellite paths at frequencies higher than 
10 GHz, which are significantly different from those at lower 
frequencies. 

Some understanding of the propagation characteristics 
has already been obtained from theoretical studies of microwaye 
transmissions at centimeter and millimeter wavelengths. A rep- 
resentative list of pertinent publications is given in Refer- 
ences 1—6. A more detailed understanding of the subject requires 
analysis and interpretation of data collected from a variety of 
experiments using satellites, rain gauges, radars, radiometers, 
terrestrial links, etc. Of all possible experiments, those using 
direct measurements of signals from a satellite aye the most de- 
sirable, since they provide a direct simulation of a communica- 
tions link along a slant path between a satellite and an earth 
terminal. 

The first satellite experiment for the determination of 
slant path signal fading statistics for frequencies over 10 GHz 
was the NASA ATS-5 15 . 3/31. 65-GHz Millimeter Wave Experiment 
(MWE) . Effects of precipitation, primarily rain, were studied 
extensively by NASA and other research organizations working with 
NASA, 7 including COMSAT Labs. The results clearly suggested that 
propagation impairments at 15.3/31.65 GHz were different from 
those at lower frequencies, and that the models established for 
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engineering applications at lower frequencies were no longer 
applicable. 

The ATS-6 20/30-GHz MWE program, developed and imple- 
mented by NASA Goddard Space Flight Center, was a follow-up and a 
refinement of the earlier ATS-5 MWE program. It provided direct 
measurements oh 20- and 30-GHz earth- satellite links from the 
eastern U.S. to the synchronous satellite located at 94 °W longi- 
tude. The satellite was launched on May 30, 1974, and became 
operational for the MWE in July 1974. It provided three modes of 
operation: continuous wave, multitone, and communications, all 

selectable at 20 and/or 30 GHz. The MWE system is described in 
detail in Reference 8. Under NASA contract NAS5-20740, COMSAT 
Labs participated in the MWE program, performing measurements of 
the 20/30-GHz satellite beacon signals and other auxiliary mea- 
surements such as radiometric sky temperature and precipitation 
at Clarksburg, Maryland. These measurements were designed to 
broaden the data base required to advance the understanding of 
propagation characteristics at frequencies over 10 GHz. 

There are several practical problems associated with a 
propagation experiment using a satellite. Due to the tremendous 
cost and complexity involved in designing and launching the 
spacecraft, assessing and assigning the priorities of various ex- 
periments operating at locations around the U.S., and scheduling 
and operating the system in coordination with various ground ter- 
minals, such a satellite propagation experiment is at most inter- 
mittent in time and spot-beam coverage, and can simultaneously 
accommodate only a few research institutions. The practical 
approach for researchers is therefore to study the results gen- 
erated by a limited number of direct satellite transmission mea- 
surements so that these results can be correlated with those 
generated by various auxiliary ground-based experiments, such as 
rain gauges, radars, radiometers, and terrestrial links. The 
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ultimata goal is to establish definitive criteria for using the 
results of the relatively inexpensive and continuously available 
ground-based experiments to make valid engineering estimates of 
the important propagation parameters for space communications 
applications . 

Such an approach is particularly pertinent to the ATS- 6 
MWE program. The MWE program was not assigned a high priority 
throughout most of the experimental period, and transmission time 
was not plentiful. In addition, the regular assigned transmis- 
sions, which averaged about one hour per day, were not always 
useful, since precipitation rarely occurred within the scheduled 
transmission time. While MWE experimenters were encouraged to 
send requests for 20/30-GHz transmission to NASA ATS-6 Operation 
Center when a precipitation event was occurring or anticipated at 
the site, responses to these requests were seldom positive due to 
other higher priority work in progress or scheduled. Obviously, 
the ATS-6 MWE was not intended to provide continuous transmis- 
sions for the determination of long-term fading statistics as 
some of the experimenters might have hoped, but rather to provide 
adequate transmissions for experimenters to study the propagation 
characteristics under various meteorological conditions based on 
a number of individual but significant events. 

In view of these factors, the COMSAT Labs MWE program 
consisted of the following measurements: satellite beacon sig- 

nals at 20 and 30 GHz; radiometric sky temperatures at 11.6, 20, 
and 30 GHz; precipitation at six sites aligned beneath the slant 
path; and rain cloud motion. Understanding the correlations 
between the beacon signals and the data from the ground-based 
measurements became a major goal of the program. 

The COMSAT experimental signal receiving configura- 
tions, including the 3- and 5-meter antenna-receiver systems, are 
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described in Section 2. Data analyses of radiometric measure- 
ments of sky temperature at 11 . 6 , 20, and 30 GHz are discussed in 
Section 3, which also includes some of the earlier sky tempera- 
ture results collected by COMSAT under the ATS- 5 MWE program for 
comparison. Section 4 is devoted to modeling of the data from 
rain gauges. A new method for using the rain gauge records to 
predict the slant path attenuation is postulated and examined. 

The satellite beacon experiment, described in Section 5, includes 
measurements of the 20/30-GHz MWE beacon at 42° elevation angle 
during the regular experimental period as well as at low eleva- 
tion angles during the period when the satellite was moving from 
its 94 °W longitude position to the eastern position at 35 °E 
longitude at the end of May 1975. Conclusions are given in 
Section 6. 

An interim contract report, ’'Preliminary 'Report on 
Atmospheric Attenuation Studies on ATS-6 Satellite 20/30 GHz 
Beacon Signals," which largely documented the hardware configura- 
tions and data analysis routines, is attached as Appendix A. Most 
of the contents of the interim report have not been duplicated in 
the present report, which focuses on the results of data analysis 
and engineering implications. 
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2 . 


SYSTEM CONFIGURATION 


Measurements were performed at COMSAT Labs, Clarksburg, 
Maryland, about 65 km northwest of Washington, D.C. The measure- 
ment systems can be summarized as follows: 

a. a 3-meter-diameter transportable antenna-receiver sys- 
tem including a 20-GIIz receiver for receiving satellite beacon 
signals and a radiometric receiving system at 11.6 GHz for ob- 
serving sky temperature using the same antenna, 

b. a 5-meter-diameter antenna-receiver system including 
20- and 30-GHz receivers for the satellite beacon signals and two 
radiometric receivers at 20 and 30 GHz for observing sky tempera- 
ture using the same antenna, and 

c. a network sf six rain gauges extending 7 km underneath 
the boresight path toward the satellite at 94 °W longitude for ex- 
amining the precipitation along the propagation path. 

In addition, a 5.5° beamwidth RCA AVQ-10 weather radar at 4.5 GHz 
was employed for large-scale visual observation of the motion 
pattern of rain clouds and to predict the approach of a precipi- 
tation cloud. 

Figures 2-1 and 2-2 are block diagrams of the 3- and 
5-meter systems, respectively. The antennas are solid-surface 
aluminum paraboloids equipped with focal-point feeds and polar- 
ization rotators. Both antennas have servo-controlled elevation 
and azimuth positioners and decimal position readouts with a 
resolution of 0.1° and 0.01° for the 3- and 5-meter antennas, 
respectively. The nominal pointing direction to ATS-6 is SSW 
204° from true north at an elevation angle of 42°. 

As indicated in the figures, tunnel diode and uncooled 
parametric amplifiers were used in the 20-GHz front ends of the 
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3- and 5-meter systems, respectively. The 30-GHz system had a 
low-noise mixer front end. IF processors which detected the 
down-converted signals by means of individual narrowband phase- 
locked receivers produced output amplitudes that were propor- 
tional to the satellite beacon signal attenuation levels. These 
outputs were then fed to stripchart recorders. Important perfor- 
mance parameters of the 3- and 5-meter systems are listed in 
Tables 2-1 and 2-2, respectively. 


Table 2-1. Important Performance Parameters of 
the 3-Meter System 


Parameter 

11.6 GHz 

20 GHz 

Antenna Gain 

48.7 dB 

53.5 dB 

Antenna Beamwidth Cl dB) 

0.3° 

0.17° 

Pointing Accuracy 

±0.1° 

±0.1° 

Noise Figure 

5 dB 

8 dB 

Radiometer Sensitivity 

1 K 

— 

Radiometer Bandwidth 

100 MHz 

«■« 

Receiver Bandwidth 

— 

10 kHz 


Table 2-2. Important Performance Parameters of 
the 5-Meter System 


Parameter 

20 GHz 

30 GHz 

Antenna Gain 
Antenna Beamwidth (1 dB) 
Pointing Accuracy 
Noise Figure 
Radiometer Sensitivity 
Radiometer Bandwidth 
Receiver Bandwidth 

57 dB 
0.12° 
±0.01° 
4 dB 
1 K 

100 MHZ 
10 kHz 

60.5 dB 
<0.10° 
±0.01° 

8 dB 
1 K 

100 MHz 
10 kHz 
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Both antenna systems were equipped with on-line radio- 
metric subsystems. Operating in a Dicke-switched mode, these 
radiometers also included a noise source at a known temperature 
for calibration purposes. The A.I.L. (Airborne Instrument Labo- 
ratories) model 239 2C radiometric receiver produced an output 
voltage proportional to the sky noise temperature. In the 
switched mode, the output was the net change of sky noise tem- 
perature from that on a clear sky day. Performance parameters of 
the A.I.L. radiometric receivers are summarized in Table 2-3, 


Table 2-3. Performance Parameters of 
A.I.L. 2392C Radiometer 


IF Bandwidth 

100 MHz (5-105 MHz) 


Noise Figure 

7 dB 


Input Impedance 

50 


Input SWR 

Less than 2:1 over the passband 

IF Gain 

85 dB maximum 


Switch Rate 

5—500 Hz, with ±0.05% 

stability 

Switched System Output 

±5 V at 2 mA 


Integration Time 

0,1, 1, 3, 5, 10, and 

30 s 


The rain gauges were tipping bucket gauges, for which 
each tip represented 0.254 mm of rainfall accumulation. The gauges 
were aligned beneath the satellite bore sight path at distances of 
0, 0.92, 1.08, 3.23, 5.53, and 7.4 km from the beacon receiver. 
Each gauge was equipped with a mechanically driven stripchart re- 
corder with a time accuracy of better than 3 minutes per week. 

Typical records of sky temperature, attenuation, and 
precipitation are shown in Figures 2-3a and 2-3b. The records of 
the 11. 6-GHz radiometer and the rain gauge at the site (0 km) 
were uninterrupted for more than a year beginning in September 
1974. Other data were also available continuously once the 
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Figure 2- 3a. Typical Sky Temperature-Attenuation Record 
from the 3-Meter System with Simultaneous Precipita- 
tion Records Taken from Ground Rain Gauges Beneath 
the Boresight Direction 0, 1.08, and 3.23 km 
Away from the 3-Meter System (elevation 
angle of antenna pointing is 42 °) 
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Figure 2-3b. Typical Records of Simultaneously Measured 
20/30-GH2 Radiometric Sky Temperatures and ATS-6 
Beacon Signals from the 5-Meter System (eleva- 
tion angle of antenna pointing is 15°) 







equipment became operational. The only exceptions were the 20/30- 
GHz satellite beacon signal data, which were rather intermittent 
due to the limited availability of the ATS-6 MWE transmission. 
Details are summarized in Table 2-4. 

In a preliminary effort, all the stripchart records, 
such as those shown in Figure 2-3, were scanned. In the subse- 
quent analysis, useless clear sky data and defective data were 
deleted first. The magnitudes were then properly scaled and the 
simultaneity in time was calculated. Next, minute-by-minute 
digitization was applied to the data before they were key punched 
on IBM cards for analysis. 


Table 2-4. Data Collection Statistics ( 


System 

Equipment 

Date Available 

3 Meter 

20-GHz CW Receiver 

July 1974 

11.7-GHz Radiometer 

July 1974 


Rain Gauge 

July 1974 


20-GHz CW Receiver 

March 1975 

20-GHz Radiometer 

March 1975 

5 Meter 

30- GHz CW Receiver 

April 1975 

30-GHz Radiometer 

April 1975 

Others 

5.4-GHz Weather Radar 

September 1974 

Rain Gauge Network 

February 1975 


ust 1974— August 1975) 

Total Data Collection Time 


150 hr 

Continuous since July 1974 
Continuous since July 1974 
30 hr 

Continuous since March 1975 
25 hr 



Continuous since April 1975 
No data collection 
Continuous since February 1975 






































3. RADIOMETRIC EXPERIMENT 


A radiometer measures the radiative thermal emission 
collected by its antenna reflector. For engineering purposes the 
emission measured by a radiometer can be separated into three 
separate terms : 

a. the transmission term (the emission from an active 
radiator, such as the sun), 

b. the absorption term (the reradiation from an atmo- 
spheric medium when it is illuminated by an active radiator) , and 

c. the scattering term (the scattering of incident elec- 
tromagnetic waves of radiation from scatterers in the medium) . 

A more detailed discussion of these terms and their precise math- 
ematical representations can be found in References 3 and 9. 

All three terms are frequency dependent. For frequen- 
cies less than 10 GHz, the scattering term is negligible since 
the single scattering albedo is less than 10 percent. 3 The ab- 
sorption term, which is absent on clear sky days, is critically 
sensitive to meteorological conditions such as clouds, rain, and 
snow. The radiation attenuation, A, which is the quantity re- 
searchers attempt to evaluate on the basis of the radiometric ex- 
periments, is included in both the transmission term and the 
absorption term. The Dicke-switched radiometer system calibrates 
an internal noise load on a clear sky day so that on any non- 
clear days the system output, compared by switching to the pre- 
calibrated internal noise load, represents the temperature 
increase from clear sky conditions contributed solely by the ab- 
sorption term. With some simplifications of the radiative heat 
transfer model, such as the assumption that the absorption is 
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uniform along a well-defined ray path through the medium, the 
following engineering equation can be derived: 

A ■ - 10 ( x - fef ) (3 - 1) 

where A = attenuation, dB 

T m eas = measured radiometric temperature, K 
Tamb = ambient temperature of the medium, K 
= 1.12 [ground temperature] - 50 

Equation (.3-1) is used to reduce the 11.6-GHz radiom- 
eter data for determining the long-term fading statistics. The 
result is shown in Figure 3-1, where T am k is assumed to be 290 K. 
Also shown in Figure 3-1 are curves known as CML curves, which 
were obtained a few years ago by COMSAT 1 0 at Clarksburg using a 
radiometer at 15.3 GHz pointing at a slightly different direction 
(antenna centered at 220° azimuth from true north at 35° eleva- 
tion toward ATS-5) . 

The difference between the two sets of curves is evi- 
dent. The CML curve was derived in two steps. The attenuation, 
A, was calculated at 15.3 GHz using equation (3-1) and then lin- 
early scaled from 15.3 GHz down to 12 GHz. Since equation (3-1) 
becomes less dependable as the frequency increases and since the 
validity of such frequency scaling is uncertain, the validity of 
the CML curves is questionable. The 11.6-GHz curve, which uti- 
lizes equation (3-1) at a relatively low frequency and involves 
no scaling, is believed to be an improved version of the earlier 
CML curves. 

The inadequacies of equation (3-1) at frequencies sig- 
nificantly higher than 10 GHz can be investigated by studying 
records of simultaneous radiometric sky temperature and satellite 
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beacon signal level measurements at 20/30 GHz. Figure 3-2 con- 
sists of the 30-GHz records and their corresponding scattergrams 
for two prominent precipitation events. Figure 3-3 is a composite 
scattergram based on data collected from several events in which 
the peak precipitation was greater than or equal to 20 mm/hr. 
Figures 3-4 and 3-5 are similar records for 20 GHz. Data points 
corresponding to signal levels less than 1 dB are omitted in the 
figures . 

As can be seen from the scattergrams, unless an ex- 
tremely low value of ambient temperature, T ain k, is assumed (which 
is rather unreasonable, since T a mb is essentially frequency inde- 
pendent) , the curves based on equation (3-1) will significantly 
underestimate the attenuation at 20/30 GHz. As the frequency 
increases, the scattering term becomes so significant that it 
effectively saturates the sky temperature measurable by a radi- 
ometer, while the satellite beacon signal attenuation level in- 
creases with frequency without such saturation. It follows that 
there is no longer any physical basis for using equation (3-1) to 
model the relationship between A and T meas . To demonstrate the 
correlations between the two quantities, simple numerical fit 
curves are provided, as shown in the figures. These curves, which 
are for empirical applications, are useful only for low to medium 
values of T am b and A. A practical engineering relationship con- 
sistent with the physics of radiative transfer at higher frequen- 
cies has yet to be derived. 

For space communications engineers interested primarily 
in long-term statistics of the signal fading levels, the empiri- 
cal method of generating the long-term fading statistics from 
radiometric data is to first establish the correlation between 
the radiometric sky temperature and signal fading level by ana- 
lyzing several reliable records obtained during prominent pre- 
cipitation events, such as the one shown in Figure 2-3a. This 
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Figure 3-2a. Simultaneous Record of Radiometric 
Sky Temperature and Satellite Beacon Signal 
at 30 GHz for an Event of June 6, 1975 
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Figure 3-2c. Simultaneous Record at 30 GHz for 
an Event of May 16, 1975 
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Figure 3-3. Composite Scattergram of All the Simultaneous 
Records of Radiometric Sky Temperatures vs Satellite 
Beacon Signals at 30 GHz {the best numerical 
fit curve is T 30 = 57.17 A 35 569 ) 
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Figure 3-4c. Simultaneous Record at 20 GHz for 
an Event of May 16, 1975 
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Figure 3-5. Composite Scattergram of All the Simultaneous 
Records of Radiometric Sky Temperatures vs Satellite 
Beacon Signals at 20 GHz (the best numerical 
fit curve is T„ n = 59.38 A?: 423 ) 



correlation is then used as a general guideline for scaling the 
large amount of radiometric data. This method is also consistent 
with the general approach of the ATS-6 20/30-GIIz I1WE, as stated 
in Section 1. 

The correlations between the 11.6-GIIz sky temperature 
and the 20-GHz beacon signal level for prominent individual 
events are plotted in Figure 3-6; Figure 3-7 is a composite plot 
for all the events. Again, data points corresponding to a beacon 
signal level less than 1 dB are omitted from the figures. Empiri- 
cal curves obtained by numerical fitting of the data are summa- 
rized in Table 3-1. Based on the composite data, the following 
relationship is derived: 

A 2Q (dB) = 7.418 x 10- 4 T}j.f* 15 13-2) 

where A 20 is the beacon signal attenuation at 20 GHz and Tu.e is 
the sky temperature measured by the 11.6-GIIz radiometer. The 
1-year Tn f{ data obtained at Clarksburg have been analyzed, and 
the final results are presented in Figure 3-8. Since the Tja.e 
reading could not be precise at high temperatures, the cumulative 
statistics for A 2 0 greater than 10 dB are shown as a dashed line. 
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a. Temperature and Attenuation vs Time for 
Event of October 16, 1974 

Figure 3-6. Simultaneous Records of Radiometric Sky 
Temperatures at 11.6 GHz and Satellite Beacon 
Signals at 20 GHz and Their Corresponding 
Scattergrams tsheet 1 of 10) 
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Figure 3-6. Simultaneous Records of Radiometric Sky 
Temperatures at 11.6 GHz and Satellite Beacon 
Signals at 20 GHz and Their Corresponding 
Scatter grams Csheet 2 of 10) 
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c. Temperature and Attenuation vs Time for 
Event of November 25, 1974 


Figure 3-6. Simultaneous Records of Radiometric Sky 
Temperatures at 11.6 GHz and Satellite Beacon 
Signals at 20 GHz and Their Corresponding 
Scattergrams (sheet 3 of 10) 
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e. Temperature and Attenuation vs Time for 
Event of March 19, 1975 

Figure 3-6. Simultaneous Records of Radiometric Sky 
Temperatures at 11.6 GHz and Satellite Beacon 
Signals at 20 GHz and Their Corresponding 
Scattergrams (sheet 5 of 10) 
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Figure 3-6. Simultaneous Records of Radiometric Sky 
Temperatures at 11.6 GHz and Satellite Beacon 
Signals at 20 GHz and Their Corresponding 
Scattergrams (sheet 7 of 10) 
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Event of Hay 22 , 1975 

Figure 3-6. Simultaneous Records of Radiometric Sky 
Temperatures at 11.6 GHz and Satellite Beacon 
Signals at 20 GHz and Their Corresponding 
Scattergrams (sheet 9 of 10] 
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Figure 3-7. Composite Scattergram of Radiometric 
Sky Temperature at 11.6 GHz vs Satellite Beacon 
Signal at 20 GHz (the numerical fit curve 
is Tn.e = 52.26 A® 6 5 4 9 ) 
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Table 3-1. Summary of T u 6 = i 


Event Date 

Figure No. 

Total 
Minutes 
of Data 

October 16, 1974 

3-6a , 3-6b 

18 

November 25, 1974 

3-6c , 3-6d 

12 

March 19, 1975 

3-6e , 3— 6f 

100 

April 29, 1975 

3-6g , 3-6h 

33 

May 22, 1975 

3~6i, 3-6 j 

35 

Composite 

3-7 

198 
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Figure 3-8. 20-GHz Cumulative Attenuation Statistics 
Derived from 11. 6-GHz Radiometric Data Based on 
Empirical Equation (3-2) 




4. RAIN GAUGE EXPERIMENT 


A basic problem in estimating the microwave attenuation 
over a satellite-earth propagation path based on rain gauge data 
is that, for a given storm event, the attenuation and the rain- 
fall records very often do not have a consistent detailed corre- 
lation. As shown in Figure 2- 3a, while the correlation between 
the radiometer curve and the satellite signal level is good, the 
rain gauge data do not have a consistent correlation with the 20- 
GHz curve. A particularly puzzling point is the peak registered 
on the rain gauge 3.23 km away from the receiver. 

The lack of consistency in correlation is not uncommon. 
For example, in an earlier ATS-5 15 . 3/31 . 65-GHz NASA experiment 
at Rosman, North Carolina, 7 it was found that neither the ground- 
averaged nor the height-averaged rainfall data of 10 rain gauges 
beneath the boresight path exhibited a sufficiently detailed cor- 
relation with the satellite signals. Hence, the correlation could 
be made meaningful only by literally shifting the time scale of 
the rain gauge data so that the peaks of the rainfall matched the 
peaks of the signal attenuation. This shifting scheme has been 
widely employed for determining the fading depth versrrs precipi- 
tation rate in both satellite beam signal and radiometric experi- 
ments. The theoretical basis for such a scheme has yet to be 
established. 

In the present COMSAT experiment, data with such appar- 
ent discrepancies in time and magnitude have been collected and 
analyzed. A new model has been established for adjusting the rain 
gauge data by taking into consideration the different falling 
times of raindrops of different sizes. This model can effectively 
account for the fact that the rain rate registered at one instant 
on a ground rain gauge away from the receiver is actually re- 
sponsible for wave attenuation at an earlier instant when the 
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rainfall intercepts the slant propagation path. For instance, 
correlations of the March 19 event shown in Figure 2- 3a can be 
improved markedly after the falling time adjustment is made. De- 
tails of the model have been published in Reference 11, which is 
included as Appendix B. 

Further improvement of the model is obviously desired. 
Most importantly, the wind motion, which certainly affects the 
falling processes of raindrops,"'' has been completely omitted. 
Evaporation, condensation, sublimation, accretion, and coales- 
cence are all involved to different degrees. In addition, it is 
necessary to consider the positions of the rain gauges with ref- 
erence to the height of the rain cloud and the characteristic 
size of the rain cells? both are dependent upon the rain type and 
other meteorological conditions. As a general rule, to avoid the 
discrepancy between the precipitation rate determined by atten- 
uation and the precipitation rate derived from the rain gauge 
fields, the rain gauges should be placed so that the eyaluation 
of the rain rate is interpolated to an elevation lower than the 
boundary of the rain cloud and so that the distance between rain 
gauges is smaller than the characteristic size of the rain cells. 
The model may also be significantly refined by using other drop 
size distributions such as the Joss distribution f which is more 
appropriate for thunderstorms in certain climatological regions 
than the Laws and Parsons distribution. 

Another aspect of the rain gauge data analysis is the 
cumulative precipitation statistics, T CR) . These statistics have 
been compiled and plotted as a solid line in Figure 4--1 based on 
records of the rain gauge for the 3-meter system from October 1, 
1974, to September 30, 1975. The annual accumulated rainfall 
depth for the period is 1100 mm. It is of interest to compare 
this curve with that derived from the Rice and Holmberg model, 12 
which has been widely used by systems engineers to scale 


4-2 


PERCENTAGE OF TIME RAIN RATE 2 ABSCISSA 





empirically cumulative statistics, with simple precipitation data 
from local weather bureau records as the only input. 

Based on past ; knowledge , 1 0 the model for the Clarksburg 
area yields 


T(R) 


100 

24 x 365 


0.033 e“ 0 • 



2 5 8 R 


+ 1.86 e“ 


0 3R + 0 . 2 (1 - 8 ) 

1 . 6 3R*] I 


(4-1) 


where 8, the ratio of annual thunderstorm rainfall depth to M, 
the total annual rainfall depth, has been designated as 0.25. 

The value of M estimated from tabulations of the local weather 
stations for the 1-year period beginning October 1, 1974, is 
1283 mm. (As a result of the saturation of the tipping bucket 
mechanism, the rain gauges used in the experiment would register 
a relatively lower volume of rain rate in heavy storm events. 
Hence, the higher value of "M 11283 mm as opposed to 1100 mm] was 
not unexpected.) The Rice and Holmberg curve of TCR) for M 
= 1283, shown as a dashed line in Figure 4-1, is reasonably close 
to the solid curve for cumulative statistics less than 0,1 per- 
cent. Incidentally, if M - 1100 in equation (4-1), the curve 
will be almost a direct overlay of the solid curve. For all prac- 
tical purposes, the Rice and Holmberg model is sufficient for 
scaling the cumulative statistics. 

It should be noted that the model 1 s annual rainfall 
depth factor, M, can have a year-to-year variation of over 
20 percent. For example, the total rainfall depth of September 
1975, a wet month, was about 250 mm more than that of September 
1974, a dry month. If the annual period is shifted one month 
forward, i.e., from September 1, 1974, to August 31 f 1975, the 
total rainfall depth will not be 1283 mm, but will be only about 
1040 mm. It therefore appears that the Rice and Holmberg model 
may not be as consistent as desired from one year to another. 
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5. SATELLITE BEACON EXPERIMENT 


As stated in Section 1, an important limitation of the 
satellite beacon experiment was the infrequency of beacon trans- 
missions during precipitation events. Although regular weekly 
transmission times averaging about one hour per day were assigned 
to COMSAT for data collection and equipment calibration, precipi- 
tation rarely occurred on-site at these times. 'Requests to NASA 
ATS Operations Control Center for 20/30-GHz beacon transmission 
when precipitation was occurring or anticipated seldom resulted 
in positive responses. There were also occasions when beacon 
transmissions were abruptly terminated while valuable experimen- 
tal data in thunderstorm conditions were being collected. 

Figure 5-1 consists of scattergrams for prominent 
events of beacon signal fading at 20 and 30 GHz. Least-mean- 
square regression fittings of the form 

A 3Q = ctA§ 0 (5-1) 

were performed. The values of a and 0 range from 1.9 to 2.18 and 
from 0.997 to 1.0.6, respectively, with the values of fitting over 
all the samples given as 

a = 1,989, 0 = 1.06 (.5-2) 

as shown in Figure 5- Id. All of these values are remarkably 
consistent with those generated by NASA based on its measurement 
of ATS-6 MWE 20/30-GHz signals at Rosman, North Carolina, and 
Greenbelt, Maryland. 8 

The equations given by Medhurst,- 13 which are improved 
versions of the equations based on Rayleigh assumptions,- 1 can. be 
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Figure 5-la. Scattergram of 30-GHz vs 20-GHz 
Attenuation Recorded May 16, 1975 
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Figure 5-lb. Scattergram of 30-GHz vs 20-GHz 
Attenuation Recorded June 7, 1975 
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used to postulate 0 the theoretical attenuation relationships at 
20 and 30 GHz as 

Azo = 0.0687R 1 - 1004 L (5-3a) 

A 3 o = 0.1649R 1 * 0 3 5 3 L (5-3b) 

where L is the equivalent path length of the precipitation 
medium. Combining equations (5-3) and (5-1) makes it possible 
to establish the empirical relationship between L and R: 

— 1 1 / (a - 3 ) 

.1 0 0 4x3-1 . 0 3 5 3 (5-4) 

for first-order estimations in engineering applications. How- 
ever, it should be noted that the concept of the effective path 
length is at most a vague one. Even in a statistical sense, the 
dependence of L on frequency, elevation angle, and other meteoro- 
logical factors such as cell size has yet to be understood. Fur- 
thermore, the’ validity of the Rayleigh assumptions at 20 and 
30 GHz is not well justified. 

In late 24ay 1975 ATS-6 began to move slowly from its 
original location toward its new position over the Indian Ocean. 
The path vector elevation angle at the COMSAT Labs 1 receiver sta- 
tion decreased steadily to zero over a 3-week period. While the 
20- and 30-GHz beacon signals were not continuously available, 
many hours of data were collected at decreasing elevation angles. 
As ATS-6 approached the horizon (local horizon was at an eleva- 
tion angle of 1.2°), strong scintillations occurred and in addi- 
tion there was an increase in the average level of the received 
signals . 

To the first order, the scintillations are believed to 
be caused by the multipath effect of the signals which rapidly 


r* 

0 


a x 0-. 0687 ^ 


0.1649 
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drift in and out of the receiving antenna beam. Since the beam- 
width of an antenna at 30 GHz is much smaller than that at 
20 GHz, the scintillations at 30 GHz are more pronounced both in 
magnitude and frequency then those at 20 GHz. This can be seen 
from Figure 5-2, which is a typical record of the simultaneous 
20- and 30-GHz scintillation measurements. 

Although antenna tracking problems and calibration 
errors, for example, might also have contributed to the signal 
variations, a continuous effort was made to minimize this effect. 
The magnitude of the scintillations showed a strong dependence on 
path elevation angle and local climatological conditions. The 
peak-to-peak magnitude of the variations increased to more than 
10 dB during the final hours of data acquisition near 3.5° eleva- 
tion angle. In addition, sharp nulls in the signal level, be- 
lieved to be caused by angle of arrival or low-elevation-angle 
multipath, were frequently observed. These scintillations had a 
periodicity ranging from 15 seconds to more than 1 minute and 
showed a large dependence on the local weather. All of these 
factors are evident from the records shown in Figures 5-2 and 
5-3. 

Based on the dynamic range of the receiver and the 
beamwidth of the antenna, peak-to-peak scintillations of 20 dB 
or more would be anticipated at low elevation angles. However, 
scintillations of this magnitude were not observed. The lowest 
angle at which the satellite made transmissions received at 
Clarksburg was 3.5°. The path was essentially a treetop path 
over a flat terrain. During the 2-day period beginning on June 
12, 1975, when the satellite descended rapidly from an elevation 
of 10° to 3.5°, there was persistent precipitation at Clarksburg 
except during the final two hours. Thus, low-elevation-angle 
effects under clear sky conditions could not be measured. In 
fact, the precipitation had contaminated the signal so that the 
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Figure 5-2. Recordings of 20/30-GHz Low-Elevation-Angle Signals 
(A and B show steady variations at 20 and 30 GHz, respectivley) 











peak-to-peak magnitudes of the fluctuations at 3.5° when the sky 
was heavily clouded but without precipitation were less than 
those at higher elevation angles (see Figure 5-4). 

The low-elevation-angle scintillations observed at 
Clarksburg are believed to be a result of the combined effects 
of angle of arrival changes, refraction and scattering by tropo- 
spheric turbulent structures, and precipitation conditions. It 
is possible to identify the effect of water-laden rain clouds 
and rain. An interesting example is provided by the 30-GHz 
recordings of Figure 5-3, which indicate that steady 5-dB peak- 
to-peak variations of signal level were continuously recorded for 
more than an hour before the scintillations gradually began to 
grow in magnitude. While the scintillation magnitude increased, 
significant absorption was also occurring, as indicated by the 
mean beacon signal strength and the radiometric output. In addi- 
tion, Figures 5-3 and 5-4 show the extreme contrast in the nature 
of scintillations recorded during rain and during early morning 
clear skies. 

Since the scintillation magnitude is expected to be 
proportional to the path length, through the lower atmosphere, the 
magnitude would be expected to be proportional to the cosecant of 
the antenna elevation angle. Because of the extreme variability 
of the local weather during the time when scintillations were 
recorded, the relationship between elevation angle and scintil- 
lation magnitude cannot be easily established. Peak-to-peak 
scintillation magnitude versus elevation angles is shown in Fig- 
ure 5-5 with the data sample taken under clear to light drizzle 
weather conditions. 

Also shown in Figure 5-5 is the ratio of the scintil- 
lation magnitudes at 20 and 30 GHz. The ratio appears to be 
relatively constant in most cases. Detailed spectrum analysis of 
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Figure 5-4. 20-GHz Recordings (a: rapid scintillations dur 
ing rainfall at ~6° elevation angle, b: strong scintilla- 
tions during heavy local cloud cover at ~6° elevation 
angle, c: slow scintillations during overcast 
skies at -3.5° elevation angle) 
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Peak-to-Peak Signal Fluctuations at 20/30 GHz 
at Low Elevation Angles 




the scintillations has not been performed due to the difficul- 
ties involved in making such an analysis from the stripchart 
data. Simple visual observation reveals that the power spectrum 
is concentrated in the 1- to 2-Hz band. Furthermore, the corre- 
lations between the 20- and 30-GHz scintillations are high. The 
normalized correlation coefficient based on evaluations of some 
short-piece data samples is about 0.8,. 
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6. CONCLUSIONS 


Unlike the earlier ATS-5 15. 3/31. 65-GIIz MWE program 
in which satellite beacon transmissions were sufficient to per- 
mit experimenters to readily establish the various cumulative 
attenuation-precipitation statistics, the ATS-6 20/30-GHz MWE 
program provided transmissions which were inadequate for directly 
deriving any long-term statistics. 

To indirectly establish the long-term cumulative 
attenuation-precipitation statistics, a detailed correlation must 
be obtained between the satellite signal data and the data gener- 
ated from auxiliary ground-based measurements on an event— by- 
event basis. The techniques for obtaining this correlation have 
been presented and employed satisfactorily in the present study. 
Furthermore, many of the indirectly derived statistics agree rea- 
sonably well (or can be reconciled) with earlier published re- 
sults generated independently. 

In particular, the COMSAT experimenters feel that the 
following improvements of the MWE program would have been ex- 
tremely beneficial: 

a. More beacon transmissions from the satellite should 
have been provided. 

b. Since digitization of the stripchart records was awk- 
ward and slow, the data processing should have been performed by 
minicomputer or microcomputer. 

c. The radar data should have been recorded for analysis. 

d. For the physics, a simple but adequate engineering 
equation similar to equation (3-1) should be derived for deducing 
satellite beacon attenuation from radiometric measurements at 
20/30 GHz. 
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e. The use of an additional receiving system with wider 
beamwidth in the low-elevation-angle measurement might have indi- 
cated whether a wide-beam antenna receiver system would improve 
satellite-earth communications on low-elevation-angle paths. 

The MWE program at COMSAT has certainly broadened the 
necessary data base and advanced the understanding of propagation 
characteristics at frequencies over 10 GHz. Additional data ob- 
tained from more refined satellite experiments such as the CTS 
and the COMSTAR programs for use . in future propagation studies 
may eliminate propagation anomalies as the major uncertainty in 
the design of communications systems at frequencies over 10 GHz. 
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APPENDIX A. PRELIMINARY REPORT ON ATMOSPHERIC 
ATTENUATION STUDIES ON ATS -6 SATELLITE 
20/30 GHz BEACON SIGNALS 



Preliminary Report On 

ATMOSPHERIC ATTENUATION STUDIES ON ATS-6 
SATELLITE 20/30 GHz BEACON SIGNALS 


D. J. Fang and J. M. Harris 
COMSAT Laboratories 
Clarksburg, Maryland 20734 

ABSTRACT 

The study is a part of a program of COMSAT praticipation in the NASA ATS-6 
millimeter wave experiment. In the initial data collection phase, the measure- 
ment equipment included: (1) A 3 meter diameter antenna-receiver system for 
accepting 20 GHz satellite CW signals and a radiometer system at 11.6 GHz 
for monitoring sky noise temperature. (2) A 4.6 meter diameter antenna-receiver 
system capable of receiving 20 and 30 GHz main CW signals plus radiometers 
at 20 and 30 GHz. (3) Six rain gauges along the bore sight path toward the 
satellite for recording instantaneous rainfall rates. (4) A weather radar at 
5.4 GHz for detecting the motion pattern of the rain cloud in large scale. 

Data analysis is done on an event-by event basis by using received satellite 
signals as basic reference for correlation with other information derived from 
the radiometers, the rain gauges and the radar. Important physical character- 
istics for prominent individual events of heavy rain storms are identified. 

I. INTRODUCTION 

There are "windows 1 ’ in the frequency spectrum above 10 GHz where water 
vapor and oxygen absorption are relatively low. Two such windows are the 
frequencies below 20 GHz and a band from about 27 to 38 GHz. It follows that 
potential for future earth-satellite communications is greater in these bands. 
However, the essential propagation data base needed by system designers for 
predicting microwave communication system margin is not yet available. 

COMSAT Labs' participation in the NASA ATS-6 Millimeter Wave Experiment 
(MWE) was aimed at collecting propagation data at Clarksburg, Maryland using 
the ATS-6 20/30 GHz beacons as the transmitting source and obtaining data for 
the long base-line diversity experiment jointly with NASA/GSFC, NRL and Westing 
house (Friendship). The data should be applicable to the general Washington, D.C. 
area and to other locations of similar climatological patterns. 
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The measurement equipment included three beacon signal receivers (two for 20 
GHz and one for 30 GHz), three radiometers (11.6, 20 and 30 GHz), six rain 
gauges beneath the boresight path toward the satellite, and one weather radar 
(5.4 GHz). The hardware aspects of these systems are discussed in Section II. 
Analysis of beacon signal data is oriented toward event-by-event analysis of 
each measured occurrence of precipitation-induced attenuation, as long term 
statistics of signal attenuation were not available because of the rarity of the 
satellite beacon transmission during precipitation events (due to the low priority 
accorded the MWE during most of the stay of ATS-6 at its western station). 

Details of the data analysis are given in Section III. 

An important facet of the work which was not anticipated in the initial study 
program was the low-elevation-angle measurement of the satellite transmission. 
This measurement was made during the period when ATS-6 moved gradually 
from its original location at 94°W Long, to its present location at 35°E Long. 

It was found that strong scientillations occurred at low elevation angles, and that 
there was a large increase of signal level as the satellite descended towards the 
horizon into the diffraction region. Preliminary analysis of the data is shown in 
Section IV. Conclusions are presented in the last section of the paper. 

II. SYSTEM CONFIGURATION 

Measurements were performed at COMSAT Labs, Clarksburg, Maryland, 
some 40 miles northwest of Washington, D.C. Figure 5-1 shows the location 
of the Clarksburg site relative to the other diversity sites at NRL, GSFC, and 
Westinghouse. The measurement system can be summarized as follows: 

(1) A3 meter diameter transportable antenna-receiver system that included 
a 20 GHz receiver for receiving satellite beacon signals and a radiometeric 
receiving system at 11.6 GHz for observing sky temperature using thd same 
antenna. 

(2) A 4,6 meter diameter antenna-receiver system that included 20 and 30 
GHz receivers for the satellite beacon signals and two radiometric receivers 

at 20 and 30 GHz for observing sky temperature using the same 4.6 meter antenna. 

(3) A network of 6 rain gauges extending for 7 kilometers underneath the 
boresight path toward the satellite for examining the precipitation under the 
propagation path. 

(4) A weather radar at 5.4 GHz for detecting the motion pattern of the rain 
cloud in large scale. 
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Figure 5-2 and 5-3 are the block diagrams of the 3 meter and 4.6 meter systems. 
The antennas are solid-surface aluminum paraboloids equipped with focal-point 
feeds and polarization rotators. Both antennas have servo-controlled elevation 
and azimuth positioners and decimal position readouts with a resolution of 0.1° 
and .01° for the 3 meter and 4.6 meter antennas, respectively. The nominal 
pointing direction to ATS-6 was 204° SSW from true north at an elevation angle 
of 42°. 

As indicated in Figures 5-2 and 5-3 tunnel diode and uncooled parametric 
amplifiers are used in the 20 GHz front-ends of the 3 meter and 4.6 meter sys- 
tems, respectively. The 30 GHz system has a low noise mixer front end. IF- 
signal processors which detect the down-converted signals by individual narrow- 
band phase-locked receivers give output amplitudes that are proportional to the 
satellite beacon signal attenuation levels. These outputs are then fed to strip- 
chart recorders or to data acquisition computers. Important performance 
parameters are shown in Table 5-1 (3m system) and Table 5-2 (4.6m system). 

Both antenna systems are equipped with on-line radiometric subsystems. Operating 
in a Dicke-switched mode, these radiometers also include a noise source at 
known temperature for calibration purposes. The A.I.L. (Airborne Instrument 
Laboratories) radiometric receivers, (Model 2392C) produces an output voltage 
proportional to the sky noise temperature. In the switched mode, the output 
is the net change of sky noise temperature from that of a clear sky day. Per- 
formance parameters of the AIL radiometric receivers are summarized in Table 
5-3. 

The rain gauges are the tipping bucket type, where each tip represents 0.01 inch 
(0.25 mm) of rainfall accumulation. The gauges are aligned beneath the satellite 
bore-sight path at distances of 0, 0.92, 1.08, 3.23, 5.53 and 7.4 kilometers from 
the beacon receiver. Each gauge is equipped with a mechanically driven strip- 
chart recorder with a time accuracy of better than 3 minutes per week. 

The radar system is an RCA AVQ-10 weather radar operating at 5.4 GHz. Three 
different ranges of 20, 50 and 150 nautical miles can be selected. Maximum 
elevation angle is 35°. A super-8mm camera can be attached to the PPI scope for 
for filming the radar echoes. Since the radar slant path is lower than the ATS-6 
slant path, no attempt has been made to incorporate the radar data for detailed 
analysis. However, the radar has proved useful as a device to predict the 
approach of a precipitating cloud. 
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Figure 5-2. Block Diagram: 20/30 GHz Receive/Radiometer System 
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Table 5-1 

Important Performance Parameters of 3m Antenna System 



11.6 GHz 

20 GHz 

Antenna Gain 

48.7 dB 

53.5 dB 

Pointing Accuracy 

±.1° 

±.1° 

Noise Figure 

5 dB 

8 dB 

Radiometer Sensitivity 

1°K 

- 

Radiometer Bandwidth 

100 MHz 

- 

Receiver Bandwidth 

— 

10 KHz 


Table 5-2 

Important Performance Parameters of 4.6m Antenna System 



20 GHz 

30 GHz 

Antenna Gain 
Pointing Accuracy 
Noise Figure 
Radiometer Sensitivity 
Radiometer Bandwidth 
Receiver Bandwidth 

57 dB 
±.01° 

4 dB 
1°K 

100 MHz 
10 KHz 

60.5 dB 
±.01° 

8 dB 
1°K 

100 MHz 
10 KHz 


Not all of the above systems were ready for measurement at the same time, as 
shown in Table 5-4. Delays in preparation of the 4.6m system were caused by 
the continuous slippage of the shipping schedule of the parametric amplifier. 

On the other hand, the systems have performed well once operational. 

III. DATA ANALYSIS 

An important cons t ra in t on the ent i re exper i men t was the lackof ava i 1 ab i 1 i ty 
of millimeter wave 20/30 GHz beacon transmissions from the satellite during 
precipitation events. The priority of the millimeter wave experiment was not 
high. There were regular weekly transmission times assigned at an average of 
about one hour per day. These transmission times, although used by COMSAT 
for data collection and equipment calibration, were not really useful due to the 
fact that precipitation rarely occurred on site at these times. Requests were 
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Table 

Performance Parameters of 

IF bandwidth 
Noise figure 
Input impedance 
Input SWR 
Gain modulation 

IF gain 

IF gain control 
Switch rate 
Switching waveform 
RF switch drive 

Switched system output 

Total power output 

Integration time 
constants: 

Switched system 

Total power system 

Monitor points 

DC reference 

Input power 
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AIL Radiometric Receivers 

100 MHz (5 to 105 MHz) 

7 dB 
50 ohms 

<2 to 1 over entire passband 

-3 to +3 dB, in 0.06-dB steps and 
continuous control of 0.1 dB 

85 dB (maximum) 

0 to -41 dB, in 1-dB steps 

5 to 500 Hz, stability ±0.05 percent 

Square wave 

0 to -14 volt square wave into 
10,000 ohms 

±5 volts maximum at 2 ma 
-5 volts maximum at 2 ma 

0. 1, 1, 3, 5, 10 and 30 seconds 

1, 5, and 10 seconds 

Second detector current; syn- 
chronous detector balance 

Up to 3.0 volts (for use with total 
power output) 

47 watt nominal (capable of 
operation at 105 to 125 and 
220 volt AC, 50 to 440 Hz). 
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System 

Equipment 

Available Since 

Comment 

Transportable 

20 GHz 

July 1974 

Varian 620 computer was equipped 

System 

CW-Receiver 


for automatic tracking and 
data collection. But it over 


11.6 GHz 

July 1974 

flowed frequently. Most time 


Radiometer 

' .. ' ' I 


only manual operation was used. 


Rain Gauge 

July 1975 


15 -Foot 
System 

20 GHz 
CW-Receiver 

March 1975 

j 

Due to interference problem 
the 20 GHz radiometer was 


20 GHz 
Radiometer 

March 1975 

. ! 

turned off when CW-Receiver 
was taking data. 


30 GHz 

April 1974 

Paramp not available 


CW-Receiver 


dynamic range 15-18 dB. 

Others 

5.4 GHz 
Radar 

September 1974 

For visual observation only 


Rain Gauges 

February 1975 

Total of 6, counting the 
one on site 













made consistently to NASA ATS Operations Control Center for 20/30 GHz beacon 
transmission whenever precipitation was occurring or anticipated. Positive 
responses to the requests were infrequent. There were several occassions 
when beacon transmissions were abruptly terminated while valuable experimental 
data in a thunderstorm condition was being collected. 

With the above constraint in mind, data analysis was performed largely on an 
event-by-event basis. This includes: 

(i) Correlation of the satellite 20 GHz beacon data and the radiometric 
sky temperature data. 

(ii) Correlation of the 20 GHz and 30 GHz signal attenuation level. 

(iii) Correlation of the satellite 20 GHz beacon signals and data from the 
rain-gauge network. 

Another phase of data collection, not included with the original plan, involved 
observation of the satellite descending toward the horizon in late May and early 
June of 1975. The low-elevation angle data will be presented in a separate 
section. 

(i) Correlation of the satellite 20 GHz beacon data and the radiometric sky 
temperature data. 

A typical data sample is shown in Figure 5-4, (the e vent of March 19, 1975, lasting 
150 minutes). The scattergram of the data is shown in Figure 5-5. Empirical 
fit curves are shown as dashed lines. Theoretically derived attenuation curves 
using linear scaling for ambient sky temperature of 278°K and 300°K are also 
shown. 

(ii) Correlation of the 20 GHz and 30 GHz signal attenuation level. 

A typical data sample is shown in Figure 5-6, the event of June 6, 1975. A 
scattergram of the data is shown in Figure 5-7. 

(iii) Correlation of the satellite 20 GHz beacon signals and data from the 
rain-gauge network. 


b 
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Figure 5-4. Typical Precipitation - attenuation reoord measured at Clarks- 
burg, Maryland. Ground rain gauges were 0 km, 1.1 km and 3.23 km from 20 
GHz receiver and 11.6 GHz radiometer. 
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20 GHz Attenuation 




Figure 5-6. Typical records of simultaneously measured 20/30 GHz 
radiometric sky temperatures and 20/30 GHz ATS-6 signals. 
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30 GHz 



20 GHz- aB 


Figure 5-7. Scattergram plot of 30 GHz Attenuation vs. 
20 GHz attenuation measured June 6, 1975. 
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A basic problem in estimating the microwave attenuations over a satellite — 
earth propagation path based on rain-gauge data is that for a given storm event, 
the attenuation and the rain-fall records very often do not have a consistent 
detailed correlation. As shown in Figure 5-4, it is seen that while the radiometer 
curve correlates well with the satellite signal level, the data from rain gauges 
do not have a consistent correlation with the 20 GHz curve. A particularly puzzling 
point is the peak registered on the rain gauge 3.23 km away from the receiver. 

We have collected and analyzed the data containing the apparent discrepancies in 
time and magnitude as those shown in Figure 5-4. A model for adjusting the 
rain gauge data is provided as follows: 

Let the rain-rate measured at an instant t by a ground rain gauge by R (t) which 
consists of raindrops of sizes from 0.05 to 0.7 centimeter with falling velocity 
V ranging from 2 to 9 meters/sec. It follows that 


R 0 (t) = Vt) (1) 

v 

where R v (t) is the portion of rain-rate falling at terminal velocity v. If the rain 
gauge is located h meters vertically below the slant propagation path, then R v (t) 
penetrates across the path at an advanced time h/v, when there is no wind. 
Therefore, the actual rain fall R h (t) at height h and instant t' is 


Rh(t') = J2 RyCO t -7 = (2) 

t. V 

The summation is over all the components of R y (t) in R c (t) for an all t and v such 
that t - h/v = t’. The decomposition of R„ into R„ for a different advanced 
time with h = 1000 meters is shown in Table 4-5. 
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Table 4-5 

Actual Rain Rate at 1000-Meter Height Estimated from a Ground Rain Gauge Rainfall Reading 
(The estimate is based on Laws and Parsons drop size distribution and Davies' terminal 
velocity profile [Medhurst, 1965 IEEE paper] 


Actual Rain Rate at 1000-m Height 
Advanced Time (minutes) 

Ground Measured 
Rain Rate (mm /hr) 

8.09 

4.14 

3.09 

2.57 

2.25 

2.07 

1.95 

1.89 

1.85 

1.83 


0.07 

0.13 

0.05 

0 

0 

0 

0 

0 

0 

0 

0.25 

0.14 

0.46 

0.39 

0.17 

0.16 

0.02 

0.01 

0 

0 

0 

1.25 

0.18 

0.07 

0.82 

0.48 

0.20 

0.08 

0.04 

0 

0 

0 

2.5 

0.24 

1.02 

1.55 

1.11 

0.59 

0.29 

0.13 

0.05 

0.02 

0 

5 

0.33 

1.44 

3.06 

3.18 

2.16 

1.26 

0.54 

0.29 

0.15 

0.09 

12.5 

0.43 

1.9 

4.6 

5.98 

4.98 

3.20 

2.05 

0.88 

0.53 

0.48 

25 

0.6 

2.7 

6.25 

9.95 

10.45 

7.8 

5.45 

3.35 

1.65 

1.8 

50 

1.0 

4.6 

8.8 

13.9 

17.1 

18.4 

15.0 

9.0 

5.8 

6.4 

100 

1.5 

6.15 

11.4 

17.55 

20.85 

26.55 

24.15 

17.85 

11.55 

12.45 

150 


L 


I 

\ 
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The model has been applied to the rain-rate records collected from the 6 rain 
gauge networks in correlating rain data with the ATS-6 20 GHz beacon signals. 

In most cases, the method has the merit of relocating and dispersing the peaks 
of precipitation curves of field rain gauges in the direction of better alignment 
with the pattern of attenuation curves. Figure 5-8 is the outcome of the adjusted 
rain-rate distribution of the same event shown in Figure 6-4. Assuming that 
the same rain rate of a rain gauge applies over the distance between the middle 
points of the two adjacent rain gauges, attenuation can then be estimated from 
theoretical considerations. The estimated value and the actual 20 GHz attenua- 
tion value are compared in Figure 5-9. The scattergram, Figure 5-9(b) derived 
from the adjusted rain data of Figure 5-8 exhibits a significant improvement in 
comparison with the scattergram, Figure 5-9(a), derived directly from Figure 5-4. 

IV. LOW ELEVATION ANGLE MEASUREMENT 

In late May, 1975 ATS-6 began to move slowly from its original location 
toward its new position over the Indian Ocean. The path vector elevation angle 
at the COMSAT Labs' receivers station decreased steadily to zero over a three 
week period of time. While the 20 GHz and 30 GHz beacon signals were not 
continuously available, several hours of data were collected at varying elevation 
angles. As ATS-6 approached the horizon, (local horizon was at an elevation 
angle of 1.2°), strong scintillations occurred as well as an increase in the 
average level of the received signals. 

The scintillations appear to be caused by irregular refraction and scattering 
in the earth's lower atmosphere. Although antenna tracking difficulties and 
calibration errors, among other things, may also have contributed to the signal 
variations, a continuous effort was made to minimize this effect. The magnitude 
of the scintillations showed a strong dependence on path elevation angle and 
local climatological conditions. Peak-to-peak magnitude of the variations grew 
to greater than 10 dB during the final hours of data acquisition near 3° elevation 
angle. In addition, sharp nulls in the signal level were frequently observed and 
believed to be caused by angle of arrival or low elevation angle multipath. These 
scintillations have a periodicity in the range of 15 seconds to greater than 1 
minute that showed a large dependence on the local weather. A typical record 
of simultaneous 20 GHz and 30 GHz scintillations is given in Figure 5 -10(a) and 
5— 10{b). 
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Model Adjusted Rain Rate 




Figure 5-8. Model adjusted rainrate from original data given in Figure 4. 
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Figure 5-9(a). Scattergram of 20 GHz attenuation calculated vs. measured for 
2 hour period on March 19, 1975 using original rain data from Figure 4. 
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Figure 5-9(b). Seattergram of 20 GHz attenuation calculated vs. measured for 
2 hour period on March 19, 1975 using model adjusted rain data from Figure 8. 
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Figure 5-10. Recordings of 20/30 GHz low elevation angle signals, (a) and (b) 
show steady variations at 20 GHz and 30 GHz respectively, (c) and (d) show 
simultaneous recordings of 30 GHz ATS-6 signal and 30 GHz radiometer as light 
rain began. Antenna elevation angle 6.5°. 














While several different physical mechanisms are probably responsible for the 
scintillations recorded, the occurrence of low elevation angle scintillation was 
expected. As noted earlier, the observed effect was no doubt a result of the 
combined influence of scattering by tropospheric irregularities, angle of 
arrival variations, rain, etc. However, most notable and identifiable was the 
effect of water-laden rain clouds and rain. The weather radar operated at the 
receiving site proved particularly useful in identifying the presence and movement 
of cloud systems. An interesting example is the 30 GHz recordings shown in 
Figures 5-10(c) and 5-10(d). Steady variations in the signal level of 5 dB peak- 
to-peak had been continuously recorded for well over an hour when the scintilla- 
tion magnitude increased, significant absorption was also occurring as indicated 
by the average decrease in the signal strength and the radiometers. Figure 5-11 
shows the extreme contrast in nature of scintillation recorded during rain and 
during early morning clear skies. 

Since the scintillation magnitude is expected to be proportional to the path length 
through the lower atmosphere, the magnitude would be expected proportional to 
the cosecant of the antenna elevation angle. Because of the extreme variability 
of the local weather throughout the time when scintillations were recorded, 
the relationship between elevation angle and scintillation magnitude cannot be 
easily obtained. In spite of the strong effect of the weather however, some ele- 
vation angle dependence is evident from the data. The peak-to-peak magnitude 
of the 30 GHz scintillations were generally greater than that simultaneously 
recorded at 20 GHz. This fact was particularly evident during cloud and rain 
scintillations. 


V. CONCLUSION 

As we have mentioned, due to the rigid assignment of priorities, there was a 
considerable amount of difficulty in securing the 20/30 GHz beacon transmission 
from the satellite when heavy precipitation was in progress at the Clarksburg 
experiment site. This severely curtailed the data base on the whole experiment. 

On the other hand, the program did yield many important results. For example, 
a new and more efficient method of estimating the slant -path microwave attenua- 
tion using a network of rain gauges was established. From event -by- event analysis, 
one can determine the accuracy of using a radiometric measurement of sky 
temperature as a substitute of direct satellite-signal attenuation measurement. 

The validity of using a radiometer to derive the yearly attenuation statistics, 
will be shown in the final report of this program when a one-year radiometric 
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Figure 5-11. 20 GHz Recordings: 

(a) Rapid scintillations during moderate rainfall at ~6° elevation angle. 

(b) Strong scintillations during heavy local cloud cover at~6° elevation angle 

(c) Slow scintillations during overcast skies at ~3.5° elevation angle. 


I 






observation has been completed. Other meaningful analyses which will be gene- 
rated later from the existing data base include an estimate of the scattering 
effect in radiometric measurement obtained by comparing the data at 11.6, 20 and 
30 GHz frequencies, and the low-ele v ation propagation characteristics obtained 
by studying the records of the descending satellite. All these will be useful in 
identifying the precipitation phenomena that degrade microwaves communications 
at 20 and 30 GHz. 


APPENDIX B. A NEW METHOD OF ESTIMATING MICROWAVE 
ATTENUATION OVER A SLANT PROPAGATION PATH 
BASED ON RAIN GAUGE DATA 
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Adjustments in time and drop size spectrum are needed for data collected 
from field rain gauges to realize the actual rainfall along the slant path 
for correlating rain gauge data with measured satellite signals. 


Introduct ion 

As a participant in the AT5S-6 20/30-GJI/. Millimeter Wave 
Experiment Program [I J, COMSAT Labs has conducted, among 
other experiments, measurements of the 20-GIi/. satellite beacon 
signals in conjunction with an 11.6-CiHz radiometer and a net- 
work of six rain gauges aligned beneath the satellite boresight 
path at distances of 0, 0.92, 1.1, 3,23, 5.53, and 7.4 km from 
the beacon receiver. A typical record of a precipitation atten- 
uation event is shown in Fig. 1. It can be seen that, while the 
radiometer curve correlates well with the satellite signal level, 
the data from the rain gauges do not have a consistent correlation 
with the 20-GHz curve, A particularly puzzling point is the peak 
registered on the rain gauge 3.23 km away from the receiver. 

The lack of consistency in correlation is not uncommon. 
For example, in an earlier. ATS- 5 35.3/31.65-GHz NASA experi- 
ment at Rosman, North Carolina [2], it was found that neither 
the ground-averaged nor the height-averaged rainfall data of 10 
rain gauges beneath the boresight path could exhibit a sufficiently 
detailed correlation with the satellite signals. Hence, the correla- 
tion could only be made meaningful by literally shifting the time 
scale of the rain gauge data so that the peaks of rainfall matched 
the peaks of the signal attenuations. This shifting scheme has 
been widely employed for determining the fading depth versus 
precipitation in satellite beam signal experiments as well as in 
radiometric experiments. The theoretical basis for such a scheme 
has yet to be established. 

Data containing the apparent discrepancies in time and 
magnitude, such as those shown in Fig. 1, have been collected 
and analyzed. The following section provides a model for 


A New Method of Estimating Microwave Attenuation over a 
Slant Propagation Path Based on Rain Gauge Data 
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Abstract — A basic problem in estimating the microwave attenuation 
over a satellite -earth propagation path on the basis of rain gauge data 
is that, for a given precipitation event, the attenuation and the rainfall 
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TABLE I 

Actual Rain Rate at 1000-m Height Estimated from Ground Rain Gauge Rainfall Reading (Estimate 
is Based on Laws and Parsons Drop Size Distribution and Davies’ Terminal Velocity Profile [5]) 


'■ _ 

Advanced Time (minutes) 

Ground Measured 
Rain Rate (mm/hr) 

S . 09 

4 . 14 

3.09 

2.57 

2,25 

2.07 

1.95 

1.89 

1.85 

1.83 


0.07 

0.13 

0.05 

0 

0 

0 

) 

0 

0 

0 

0.25 

0.14 

0 . 4 G 

0.39 

0.17 

0.06 

0.02 

0.01 

0 

0 

0 

1.25 

0.18 

0.70 

0.82 

0.48 

0.20 

0.08 

0.04 

0 

0 

0 

2.5 

0.24 

1.02 

1.55 

1.11 

0.59 

0.29 

0.13 

0,05 

0.02 

0 

5 

0.33 

1.44 

3.06 

3.18 

2.16 

1.26 

0.54 

0.29 

0.15 

0.09 

12.5 

0.43 

1.9 

4.6 

5.98 

4.98 

3.20 

2.05 

0.88 

0.53 

0.48 

25 

0.6 

2.7 

6.25 

9.95 

10.45 

7.8 

5.45 

3.35 

1.65 

1.8 

50 

1.0 

4 . 6 

8.8 

13.9 

17.1 

18.4 

15.0 

9.0 

5.8 

6.4 

100 

1.5 

6.15 

11.4 

17.55 

20.85 

26.55 

24.15 

17.85 

11.55 

12.45 

150 


adjusting the rain gauge data by taking the falling time of the 
raindrops into consideration. As shown in the third section, the 
adjusted rain data will provide a remarkably improved correlation 
with the signal level. Some of the defects and possible future 
refinements of the method will be discussed in the last section. 

Rain Drop Sizes and Terminal Velocities 

The basic concept that rain drops take time to fall to the 
ground is intuitive [3], [4]. The rain rate registered at one 
instant on a ground rain gauge away from the receiver is actually 
responsible for wave attenuation at an earlier instant when the 
rain fall intercepts the slant propagation path. The detailed 
physics is not so intuitive because the drop size distribution is a 
function of rain rate and drops of different sizes fall at different 
terminal velocities. 

Discussions of various drop size distributions with reference 
to different types of precipitation and of terminal velocities 
versus drop size are lengthy and beyond the scope of this paper. 
To illustrate the formulation of the above basic concept for 
actual application, it is sufficient to adopt the widely used 
Laws and Parsons drop size distribution and Davis’ terminal 
velocity profile [5]. 

Let the rain rate measured by a ground gauge at instant t be 
R G (t ) for 0.05- to 0.7-crn raindrops with falling velocity v 
ranging from 2 to 9 m/s. It follows that 

Rcit) = I Rvit) (1) 


where R y (t) is the portion of rain rate falling at terminal velocity 
v. If that rain gauge is located h meters vertically below the slant 
propagation path, then R y (l) penetrates across the path at an 
advanced time hjv when there is no wind. Therefore, the actual 
rainfall R h (t) at height h and instant t' is 

( 2 ) 

r-Wt=i‘ 


R„(n = 

f.U 


The summation for R h at instant t' is over all the components of 
R v (t) in R a (t ) for all t and v satisfying the relation ( t — hjv) ~ 
The decomposition of R G into R y for a different advanced time 
with h = 1000 m is shown in Table I. 


MODEL adjusted rain rate 



1735 1.800 1825 

TIME (GMT) (MARCH 19. 1975) 

Fig. 2. Model adjusted rain rate intercepting slant path of 42° elevation 
angle based on original data given in Fig. 1. 


The procedure of adjusting the surface rain rate R a to obtain 
the actual, rain rate R h responsible for microwave attenuation 
can be illustrated as follows. Assume that the elevation angle 
of the slant path is 30° and a field rain gauge 1.2 km away 
from the receiver (ft = 693 m) registers 100 mm/h at instant t. 
From Table I it can be determined a 13.9-mm/h rate crossed 
the slant path at an advanced time of 2.57 x (693/1000) min. 
The 13.9-mm/h rate is responsible for microwave attenuation 
at the receiver at time t — 1.78 min, etc. In actual data proces- 
sing, the sampling rate of the rain rate should be about I sample 
per second. The net effect of the adjustment is not only an 
advancement of the time of R a , but also a dispersement of the 
rain rate peak into a less steep profile. Such an effect has been 
frequently observed when comparing the attenuation data to the 
data from a far-field rain gauge. 
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ATTENUATION MEASURED (dB) 

(c) 

Fig. 3. Scattergram of 20-GHz attenuation calculated versus measured for 
2-hour period on March 19, 1975. (a) Using original data from Fig. 1. 
(b) Using original rain data shifted so that peak rain rates occur at time 
of peak measured attenuation, (c) Using model adjusted rain data from 
Fig. 2. 


Discussion 

The method postulated here has been applied at COMSAT 
Labs to the rain rate records collected from its network of six 
rain gauges for correlation with the ATS-6 20-GHz beacon 
signals. In most cases, the method does have the merit of re- 
locating and dispersing the peaks of precipitation curves of field 
rain gauges for better alignment with the pattern of attenuation 
curves. Fig. 2 is the outcome of the adjusted rain rate distribu- 
tion of the same event shown in Fig. I. If it is assumed that the 
rain rate measured at a rain gauge applies over the distance 
between the middle points of the two adjacent rain gauges, 
attenuation can then be estimated from theoretical consider- 
ations [6], The estimated value and the actual 20-GHz attenua- 
tion value are compared in the scattergrams presented as Figs. 
3(a)— 3(c). Fig. 3(a) is derived directly from the data shown in 
Fig. 1, while to derive Fig. 3(b), the measured data have been 
literally shifted, forward or backward in time, so the that peak 
rain rates occur at the time of peak attenuation. Fig. 3(c), which 
is obtained by using the model adjusted rain rates of Fig. 2, does 
exhibit a significant improvement in correlation relative to 
Figs. 3(a) and 3(b). 

To correlate the slant microwave attenuation with rainfall 
from a field rain gauge, researchers often postulate a vertically 
decreasing rainfall profile. For the NASA-Rosman ATS-5 
experiment [2], the profile was assumed to be linear with a 
height of zero precipitation assigned at 3.22 km. It is worthwhile 
to note that the present method reveals the physics of decreasing 
rainfall magnitude due to the readjustment of raindrops of dif- 
ferent sizes, with increasing height, As can be seen from Fig, 2, 
the profile is not linear, and it depends not only on the magni- 
tude. but also on the time distribution of the surface rain rate. 

In a comparison of rain rate data derived from radars (3 and 
8 GHz) ana field rain gauges [7], it was found that there were 
distinctive time delays between them. A typical delay was about 
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4 min for a rain rate less than 10 mm/h registered by a gauge 
located about 900 m beneath the radar path and 3 min or less 
for ruin rates higher than 10 mm/h. Such time differences cun 
be easily predicted by using the present method if the overall 
advanced time averaged over the percentages of rainfall dis- 
tribution is evaluated from Table I. 

CONCLUSION , , 

The method presented here is intended to improve the corre- 
lation between the data from rain gauges and the microwave 
attenuation caused by rain. Only with good correlation is it 
possible to asses the validity of using rain data as a prediction 
of satellite signal attenuation level. Needless to say, further 
improvement of the method is needed. Most important of all, the 
wind motion, which certainly affects the raindrops’ falling pro- 
cesses, has been completely omitted. Evaporation, condensation, 
sublimation, accretion, and coalescence are all involved to 
different degrees. 

In addition, it is necessary to consider the positions of the 
rain gauges with reference to the height of the rain cloud and 
the characteristic size of the rain cells; both are dependent upon 
the rain type and other meteorological conditions. As a general 
rule, in order to avoid the discrepancy between the precipitation 
rate determined by attenuation and the precipitation rate derived 
from the rain gauge fields, the rain gauges should be placed so 


that the evaluation of R h is at a height lower than the boundary 
of the rain cloud and so that the distance between rain gauges is 
smaller than the characteristic size of the lain cells. 

All of these factors are extremely complicated, and efforts 
to incorporate them into precipitation attenuation studies have 
yet to be made. The method may also be significantly refined by 
using other distributions such as the Joss distribution [8], which 
is more appropriate for thunderstorms in certain climatological 
regions, instead of the Laws and Parsons distribution. 
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